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Abstract 

Transformation matrices between different coupling schemes are required, if a reliable 
classification of the level structure is to be obtained for open-shell atoms and ions. While, 
for instance, relativistic computations are traditionally carried out in jj— coupling, a LSJ 
coupling notation often occurs much more appropriate for classifying the valence-shell 
structure of atoms. Apart from the (known) transformation of single open shells, however, 
further demand on proper transformation coefficients has recently arose from the study of 
open d— and /—shell elements, the analysis of multiple-excited levels, or the investigation 
on inner-shell phenomena. Therefore, in order to facilitate a simple access to LS <-> jj 
transformation matrices, here we present an extension to the Racah program for the 
set-up and the transformation of symmetry-adapted functions. A flexible notation is 
introduced for defining and for manipulating open-shell configurations at different level 
of complexity which can be extended also to other coupling schemes and, hence, may help 
determine an optimum classification of atomic levels and processes in the future. 



PACS: 3.65F, 2.90+p. 
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PROGRAM SUMMARY 



Title of program: Racah 
Catalogue identifier: ADQP 



Program Summary URL: http://cpc.cs.qub.ac.uk/summaries/ADQP 



Program obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland. 

Licensing provisions: None. 

Computers for which the program is designed: 

All computers with a license of the computer algebra package Maple [1]. 
Installations: University of Kassel (Germany). 

Operating systems under which the program has been tested: Linux 6.1+. 
Program language used: Maple V, Release 6 and 7. 
Memory required to execute with typical data: 30 MB. 

Keywords: angular momentum, complex atom, configuration state function, jj— coupling, 
LS— coupling, LS — jj transformation, LSJ and jj J spectroscopic notation, nonrelativistic, 
relativistic, subshell state. 

Nature of the physical problem: 

For open-shell atoms and ions, a reliable classification of the level structure often requires the 
knowledge of the LS—jj transformation matrices in order to find the main components of the 
wave functions as well as their proper spectroscopic notation. Apart from the transformation 
of individual (sub-) shell states, matrices of much larger complexity arise for the transfor- 
mation of symmetry-adapted configuration state functions which are constructed from the 
coupling of two or more open shells. 

Method of solution: 

LS — jj transformation matrices are provided for all (sub-) shell states with orbital angular 
momenta I < 3 in the framework of the Racah program [2]. These matrices are then utilized 
to transform symmetry-adapted configuration state functions (CSF), including the coupling 
of two open shells. Moreover, a simple notation is introduced to handle such symmetry 
functions interactively and to transform even atomic states which are given as a superposition 
of CSF. 

Restrictions onto the complexity of the problem: 

The program presently supports all shell states with I < 3, i.e. up to open /—shells, in 
LS*— coupling and with j < 7/2, i.e. up to open fa/2— and 57/2— subshells, in jj— coupling. 
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For the transformation of configuration state functions, the coupling of two open LS— shells 
or, correspondingly, four jj— subshells are also supported. In jj— coupling, however, a stan- 
dard order [cf. subsection 2.2] is always assumed for the coupling sequence of the individual 
shells. Several simplifications are used on the basis of this standard order. 

Unusual features of the program: 

Apart from the "interactive access" to the LS — jj transformation matrix elements between 
(sub-) shell states in LS— and jj— coupling, a complete transformation of the coupling scheme 
can be carried out also for configuration respectively atomic state functions, just by typing 
a few lines at Maple's prompt. To simplify the handling of the program, a short but very 
powerful notation has been introduced which help the user to construct stepwise symmetry- 
adapted functions of different complexity. But although the program presently supports only 
shell states in LS— and jj— coupling, the same notation can be extended also to incorporate 
further coupling schemes in the future. The main commands of the present extension are 
described in detail in Appendix B; for a quick reference on the current capabilities of the 
Racah program, we refer the reader to Ref. [3] and to a list of all available commands in the 
file Racah-commands .ps which is appended to the code. 

Typical running time: The program replies promptly on most requests. Even large tabula- 
tions of LS — jj transformation matrices can be carried out in a few (tens of) seconds. 

References: 

[1] Maple is a registered trademark of Waterloo Maple Inc. 

[2] S. Fritzsche, Comp. Phys. Commun. 103, 51 (1997); G. Gaigalas, S. Fritzsche, B. Fricke, 
Comp. Phys. Commun. 135, 219 (2001). 

[3] S. Fritzsche, T. Inghoff, T. Bastug and M. Tomaselli, Comp. Phys. Commun. 139, 314 
(2001). 



3 



LONG WRITE-UP 



1 Introduction 

The classification of the level structure of open-shell atoms and ions is a non-trivial task 
which occurs frequently in the interpretation of complex spectra. In the analysis of optical 
spectra, for instance, the correct knowledge of the LSJ spectroscopic notation of the atomic 
states may help isolate individual levels and terms without that the theoretical energies from 
ab-initio computations need be accurate enough for a direct assignment of the observed lines. 
In fact, such a demand arises already for rather simple shell structures such as the spectrum 
of Ne II pQ, for which the lowest excited 2s2p G 2 S term occurs high-up in the theoretical 
level structure, even if a sizeable wave function expansion is applied, and therefore may lead 
to misassignments — if no additional information about further properties of these levels 
or about their representation in different coupling schemes is available. Since, today, most 
relativistic computations are carried out in jj— coupling, an efficient and reliable LS — jj 
transformation of atomic states is of primary interest. 

For atoms with a single open shell and, in particular with an open s— or p— shell, the LS — jj 
transformation matrices are well known and can be obtained from different sources [2J 11] . 
These matrices are also the building blocks for the transformation of all symmetry-adapted 
functions and are often simply abbreviated by 



i.e. in terms of Fourier coefficients of the corresponding shell states in the (re-coupled) basis. 



However, such matrices for a single-shell configuration with fixed occupation N = N + N 
are only of little help in transforming atomic or configuration states with a more complex shell 
structures for which the individual (sub-) shell states need to be treated consistently with 
respect to their definition, choice of quantum numbers as well as their phase relation to each 
other. Therefore, in order to extent the single-shell matrices Q to open d— and /—shell con- 
figurations or to evaluate these transformation matrices for complex shell structures, insight 
into the construction of the subshell states is required for the all the coupling schemes under 
consideration. Moreover, the size of the transformation matrices (^Q) increases rapidly with 
the orbital angular momentum I due to the large number of (allowed) projections mi of the 
electrons in any open-shell configuration {n\li) Nl , (JI2I2) , ••• 0- For these two reasons 
and due to the complexity of the recoupling coefficients, which arise in the evaluation of the 
transformation coefficients [cf. Eqs. (|22H23|) ]. these matrix elements are often not available 
from the literature, not to mention the efficiency of their use if more than one open shell is 
involved and if such transformation need to be carried out explicitly. 

Today, an alternative and much simpler access to the transformation between different cou- 
pling schemes is possible by means of computer-algebraic manipulations. For the coupling of 





+ 
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angular momenta, for instance, such a framework for symbolic manipulations have been de- 
veloped by us during recent years and is now known as the Racah program 6 . This program 
is a powerful tool in simplifying formal expressions from the theory of angular momentum. 
Recent developments to this package concerned not only the fast and reliable evaluation of 
Racah expression but also the implementation of standard quantities , spherical harmonics 
[S] as well as the evaluation of recoupling coefficients jHj- Therefore, the Racah package also 
meets (most of) the basic requirements which are needed for the transformation of general, 
symmetry-adapted functions between different coupling schemes. With the present extension 
to the Racah package, we now support a convenient set-up and application of the LS — jj 
transformation matrices for all atomic (sub-)shells with orbital angular momenta I < 3. In 
addition, a powerful notation is provided for dealing with symmetry-adapted functions at 
different level of complexity such as atomic and configuration state functions as obtained from 
relativistic computations. In the present implementation, we support the transformation of 
such symmetry functions with up to two open (nonrelativistic) LS— shells or up to four (rel- 
ativistic) jj— subshells, respectively. For even more complex shell structures, moreover, we 
intent to utilize and implement these developments directly into the available atomic codes 
such as Grasp92 10 or the Ratip package But already with the present extension of 
the Racah program, a major step in the LSJ classification of atomic and ionic levels has 
been achieved. 

In the next section, we first explain the construction of symmetry-adapted functions in LS— 
and j j— coupling, respectively, as well as the evaluation of the transformation matrices. This 
is followed in section 3 by a short review about Racah's program structure and how it is 
distributed before we illustrate and discuss several examples in section 4. Apart from the 
transformation of configuration states with a single open shell, our third example displays the 
transformation of two atomic levels as they may arise in standard computations. Section 5 
outlines the algebraic evaluation and simplification of the transformation matrices, of course, 
by making use again of Racah itself. This section points to the road which we need to 
go in order to deal with general open-shell states and their transformation among different 
coupling schemes in the future. Finally, a few comments on further and highly desirable 
extensions of the present work are given in section 6. 

2 LS — jj transformation matrices 
2.1 Transformation of subshell states 

For a successful transformation of symmetry-adapted functions from one coupling scheme to 
another, it is first necessary to understand the construction of these functions in some detail. 
In atomic shell theory, symmetry-adapted configuration states are usually constructed from 
antisymmetrized states of N equivalent electrons of a given shell (nl), to which we briefly 
refer as (sub-) shell states below. In LS— coupling, for example, such a subshell states of the 
shell (nl) is written as [Hj 




(2) 
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where a represents all additional quantum numbers which, apart from the total orbital angu- 
lar momentum L and total spin S, are needed for the unique classification of these states. In 
practise, an additional number a is needed only for subshells with orbital angular momenta 
I > 3, i.e. for electrons from the /-, g-, . . . shells. A list of all possible subshell states for 
open s-, p-, d- and /-shells, both in quasispin and seniority notation, has been displayed 
previously in Ref. JZj, table 1. For the subshell states (J2J), moreover, the angular momenta L 
and S can be coupled also to an total angular momentum J , nl N aLSJ^j , which gives rise 
to the so-called LSJ notation. Of course, further additional intermediate angular momenta 
will arise if the subshell states of two or more open shells are coupled to each other which, 
however, should not be confused with the current discussion about the subshell states for 
equivalent electrons. 

In jj— coupling, similarly, the subshell states of N equivalent electrons of a subshell {nn) are 
represented by 

nn N uj) (3) 



where k is the relativistic (angular momentum) quantum number 

K = ±(j±l/2) for I = j ±1/2 (4) 

and two further quantum numbers v and J are found sufficient to classify all subshell states 
with j = 1/2, 3/2, 5/2, and 7/2 unambiguously. In this coupling notation, an additional 
quantum number, a, only occurs for subshell states with j > 9/2; for j = 9/2, we use the 
quantum number w = 0, 1, or 2 similar as for /—shells in LS— coupling. A list of all allowed 
subshell states in jj— coupling with j = 1/2, 3/2, 5/2, 7/2, and 9/2 were given in Ref. 0, 
table 2. In fact, all LS jj transformations of symmetry-adapted functions can always 
be traced back to the corresponding transformation of the subshell states and (j3J), from 
which these symmetry functions are built-up. 

Although, at the first glance, the definition of the subshell states in LS— and jj— coupling 
appears very similar, these states generally belong to different irreducible representations of 
the SO3 rotation group. In jj— coupling, each (nonrelativistic) reZ-shell is usually 'separable' 
into two (relativistic) subshells with total angular momenta j± = Z ± 1/2. Therefore, in order 
to transform a shell state I N aLS) into a jj— coupled basis 1 , two subshell states with j_ and 



j + may both occur in the expansion, i.e. 



-N-\ 

K u,J ) and 



+ 
N 



vj), where we utilize again the 



relativistic quantum number k to simplify the notation below. Obviously, also, N = N + N 
and k= — (k +1) > must hold where the notation k > and k < becomes clearer 
when one considers the corresponding total angular momentum j± = I ±1/2. 



1 Here and in the following, we often omit the principal quantum number n in the notation of the subshell 
states as this quantum number is irrelevant for the transformation properties of these states. The principal 
quantum number is needed only if two or more subshell states with the same I and j but different n's later 
occur in the construction of the symmetry functions. 
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Making use of this notation, the transformation between the subshell states in LS— and 
jj— coupling can be written as 



\l N aLSJ) = 



---++ 

Nv JuJ 



_N_- + (N-N) + + \ / _N_- + (N-N) + + \ 

(k vJ, k uJ) J ) ( (k vJ, k vJ) J\ I aLS J ) (5) 



+ 

_N_- ,N 



(fC V J, ~K tj) J j = 



*LS 



: ( ™ aLS j) I 1 ^ aLS J \ (fvJ, f $J) j\ (6) 



which, in both cases, includes a summation over all the quantum numbers (except of k and 
I). Here 

J which is built from the corresponding jj— coupled subshell states with j± = I ±1/2 and 



_N_- + N + + \ 

k vJ, k vJ) J ) is a coupled state with well-defined total angular momentum 



the total subshell angular momenta J and J, respectively. 
An explicit expression for the transformation coefficients 



+ 



/ _N_- j_(N-N),+ \ I ~ + _N - ,N,+ \ 

Uk uJ,k tj) J\l N aLS Jj = (l (N+N) aLS J\(k vJ, k vJ) Jj (7) 

in (JSJ) and © can be obtained only if we take the construction of the subshell states of 
N equivalent electrons from their corresponding parent states with N — 1 electrons into 
account. For a number of special configurations, expressions for these coefficients have been 
displayed before in Ref. [121 ■ I* 1 general, however, the recursive definition of the subshell 
states, out of their parent states, also leads to a recursive generation of the transformation 
matrices which we summarize in Appendix A. For the moment, it is sufficient to say 
that these transformation coefficients can be chosen real and that they occur very frequently 
as the building blocks in the transformation of all symmetry functions. The transformation 
matrices (J7| are therefore implemented in a suitable form for all (sub-) shells with / < 3 and 
occupation numbers iV = 1, 2, 21 + 1 in the current extension to the FlACAH program. 
For all other allowed occupation numbers N = 21 + 2, 4/ + 2, these transformation 
coefficients are obtained according to their electron-hole symmetry from the matrix elements 
for N' = 4/ + 2 — N. Such a symmetry relation was established originally by Grant et al. 
[PS] and later utilized also in the tabulations of Gaigalas et al. [3J for all subshells with / < 3. 

2.2 Coupling of subshell states 

Of course, many electron configurations with a single open shell occur in the notation of 
atomic levels and may allow a rough characterization. For a detailed representation of these 
levels, however, configuration state functions (CSF) with several open shell need to be taken 
into account, a situation which is even strongly enhanced when open d— or /—shell elements 
or excited levels are to be considered. In such situations, the construction of a suitable 
symmetry-adapted basis for the representation of the atomic states also requires the coupling 
of two or more (open) subshell states. In LS— coupling, typically, a CSF basis is constructed 
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from a stepwise coupling of the individual shells 1 , l 2 2 , ... 

(...(((Z^axL^i, l 2 N2 a 2 L 2 S 2 )L 12 S 12 , lf s a 3 L 3 S 3 )L 123 S 123 )...) j) (8) 

which could be written explicitly also in terms of a Clebsch-Gordan expansion. For the case 
of two open shells, for example, a proper CSF basis 



(Z x 1 aiLi5i, £ 2 2 a 2 L 2 S 2 )LS J 

^axL^M^Ms^ l^ 2 a 2 L 2 S 2 M L2 M s , 

M Ll M Sl M L2 Ms 2 M L M s 

x (L ± M Ll L 2 M L2 | LMl) (StMs, S 2 M S2 \ SM S ) (LM L SM S \ JMj) (9) 

includes the coupling of the subshell orbital angular momenta L\ and L 2 to a total L and 
the subshell spins S\ and S 2 to a total S which are finally coupled to a total J. 

.N 2 



A very similar sequence for the coupling of the subshell states 
applied also in jj— coupling. In principle, again, any (predefined) sequence of the jj— coupled 
subshells will give rise to a valid many-particle basis. For practical purposes and in particular 
for an efficient transformation of such configuration states, however, it is useful to define a 
standard order for jj— coupled configuration states such as 

+ - + + \ 

_JVi_ - + N 1+ + _N 2 _ - f +N 2+ + _JV 3 _ - f +N 3+ + \ 

(...(((((ft! U1J1, Kj I>l«7l)Jl, K 2 ^2^2)^12, «2 VlJl)J\2 ^3 ^3 ^3 ) Jl23 1 K 3 V3Ja)Jl23—)J / ( l0 ) 

which fulfills two additional conditions: 

(i) If both subshells with common /j, i.e. k- l and Kj appears in the expansion, these 

Ni_ - + N l+ + 

two subshells always occur successively in the sequence {n i ViJ%, K i f jJj) J{. Formally, 
we can use this sequence even for subshell states with zero occupation if we interprete 
\k°v = J = 0) = 1; in this case, the full Clebsch-Gordan expansion [cf. ([9*])] remains valid 
due to the orthonormality properties of the Clebsch-Gordan coefficients. 

(ii) For the LS — jj transformation of configuration states 

Inn n -Ni- - + 

( 1 a 1 L 1 S 1 , l 2 2 a 2 L 2 S 2 )Li 2 S 12 , l 3 3 a 3 L 3 S 3 )L 123 S 123 )...) J\ (...(((((k 1 v%Jx, K t v 1 J 1 )J 1 , 



_W2_ - + N 2+ + _N 3 _ - i + N 3+ + \ 

K 2 U 2 J 2 )j' l2 , K 2 V 2 3 2 )3 12 K 3 V 3 J 3 )J l23 , K 3 ^3^3)^123— ) J ) (H) 

we further assume in standard order that l\ =K\= —{k\ +1), l 2 =k 2 = —(k 2 +1), 
i.e. that the sequence of (sub-)shell states is the same on both sides of the transformation 
matrix. 

In the following, we always assume this standard order in the derivation of explicit expressions 
for the LS — jj transformation matrices. The two conditions (i — ii) simplifies the implemen- 
tation of these matrices considerably without much loss of generality. The same or at least 



S 



a very similar order in the sequence of the individual shells is assumed in many standard 
programs on atomic structure. 

We are now prepared to write down the transformation coefficients 

- Nl - - + Nl + + - N ' 2 - - +^ 2 + + N N \ 

(((«! v\J\, K x viJ\)J\, k 2 v 2 .] 2 )J 12 , n 2 ^2^2) J I d 1 aciLiSi, l 2 2 a 2 L 2 S 2 )LSJ ) (12) 

for the coupling of two open shells in LS— coupling or up to four open subshells in 
jj— coupling, respectively. As before, an expansion of a LS— coupled CSF in an appropriate 
jj— coupled basis is obtained by a summation over all quantum numbers apart from the Ki's. 
Making use of the two expansions JSJ , © and of the recoupling theory of angular momenta, 
we easily find 



(l x Nl a x L x S x , l^ 2 a 2 L 2 S 2 )LSJ 



N 2 , 



E 



JVi JV 2 v 1 J 1 v\ Ji v 2 J 2 v 2 J 2 J% J 2 J' 12 



JVi. 



(((«! V X J X , K 



(N 1 -N 1 ) + + 



N 2 



V X J X )J X , K 2 V 2 J 2 )J[ 2 , K 



(N 2 -N 2 ) + + 



V2J2) J 



x (((L X ,L 2 )L,(S X ,S 2 )S)J\ {{L X ,S X )J X ,{L 2 ,S 2 )J 2 )J) 



x ( ((J l5 J x )Ji, (J 2l J 2 ) Ja) J I (((A, A) Ji, J 2) J( 2 , J2) J 



x ( (kj ^1 Ji, /% ^1 J\)Ji I ^ 1 J x 



/ _^2_ - +(^2-^2)+ + 

x ( (k 2 !/ 2 J2, k 2 i/ 2 J2) J 2 I ' 2 2 OL 2 L 2 S 2 J 2 



(13) 



where the third and fourth line denote two recoupling coefficients which ensure that the 
sequence of the couplings on the left- and right-hand side of the expansion is taken into 
account properly. These recoupling coefficients can be evaluated by means of the Racah 
program; here we postpone the derivation until section 5 and display only the final result 

(l x Nl a x L x S x , l 2 2 a 2 L 2 S 2 )LS J) 



E 

----++--++ 

N^N 2 vi J\v 1 J xv 2 J 2 v 2 J 2 JiJ 12 



_JVi_ - + (N 1 -N 1 ) + + _N 2 _ - + (N 2 -N 2 ) + + \ 

(((K X VlJl, K X V X J X )J X , K 2 V 2 J 2 )J' X2 , K 2 V 2 J 2 )J 



Ni. 



x (-\y*+Wi+J J[J X ,J[ 2 ,L,S\ ( {k x v x J x , k x v x J x ) Ji \ l x Nl a x L x S x J x 



■h 



L x S x J x 
L 2 S 2 J2 
L S J 



Ji J2 Ju 
J 2 J J 2 



_N 2 _ - 



,N2_ 



(k 2 v 2 J 2 , k 2 v 2 J 2 )J 2 I l 2 2 a 2 L 2 S* 2 J 2 



(14) 
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Explicit representations of the transformation matrices becomes quickly cumbersome if more 
than two open shells are involved. For such complex shell structures, a symbolic and au- 
tomatic treatment seems to be inevitable. In the present implementation of the Racah 
program, an automatic transformation of configuration state functions is also restricted to 
symmetry functions with two open (LS-) shells, see section 4. 



2.3 Transformation of atomic states 

Any set of configuration state functions (JBJ) or (jlOj) can be utilized to form a many-particle 
basis for the set-up of the (corresponding) Hamiltonian matrix. The diagonalization of this 
matrix then results in approximate atomic states, given as a superposition of the accordings 
CSF. While LS— coupled CSF Q are frequently applied in nonrelativistic atomic structure 
calculations, the jj— coupled CSF (jlOj) are the basis of most relativistic computations. How- 
ever, not much need to be said here about the details of atomic structure theory. For the 
sake of simplicity, we can restrict ourselves to the configuration interaction (CI) approach in 
which (any) approximate atomic state is either written in terms of a LS— coupled 

* a (J P )) = J2 a r LS) (*)hrLS JP) (15) 



*^ P )) = E \^sJP) , (16) 



or j j— coupled basis 



where a = 1, 2, ... enumerate the atomic levels (of the given symmetry) and 7 r , ^y s denote 
the sets of all quantum numbers as required for an unique classification of the CSF. No 
further information is needed to understand the transformation behaviour when one wants 
to transform a state 



^a(J )) from one into another basis: 



Obviously, the transformation of any atomic state can be reduced simply to the transforma- 
tion of the underlying configuration state functions LS JP) and |7 S JP) , respectively. 
For these CSF, explicit expressions have been discussed in the previous subsection. 

From the user's viewpoint, of course, the main request to the present implementation concerns 
the transformation of atomic states between a LS— and jj— coupled basis. For these two 
coupling schemes, both directions of the transformation are equally supported, including 
a LS— and SL— coupling sequence. To facilitate the use of the program as well as the 
communication inside of the code, a powerful notation has been introduced to define and 
to manipulate atomic states of type (|T5|) or ((TH|) : see section 4 and Appendix B.l below. 
With this additional notations, the transformation of atomic states is achieved by a single 
line which can be typed in interactively. But already from the notation in this section it 
becomes clear how (easily) other coupling schemes such as JK— or LK— coupling ^2] could 
be incorporated into the program, for instance, in order to search for the best representation 
of some level(s) in terms of a single coupling scheme. 
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Table 1 Additional commands to the Racah package for LS — jj transformations of coupled (sub-) 
shell states. A more detailed description of these new procedures, which are provided for interactive 
work, is given in Appendix B. Our new commands also include a number of auxiliarity procedures to 
facilitate the communication with and within the Racah program. 



shell_jj(), shelLLS() 
csf_jj(), csf_LS() 

Racah_matrix_LS_jj () 



Racah_transform_csf ( ) 



Racah_transform_asf ( ) 



Racah_set .coupling ( ) 



To represent a (sub-) shell state in either jj— or LS— coupling. 

To represent a configuration state functions of one or two open shells 
in either jj— or LS"— coupling. 

Returns the LS — jj transformation matrix 



JV 



+ 
JV 



( I vLS J | (k vJ, k vJ) J J for all subshclls with / = 0, . . . , 3. 

Transforms a configuration state function from jj— to LS— coup- 
ling or vice versa. CSF of different complexity are supported by the 
program. 

Transforms an atomic state function (ASF) from jj— to LS— coupling 
or vice versa. 

Defines the coupling order in the set-up of the LS— shell states to 
either LS— or Si— coupling. 



3 New procedures for the Racah package 

Since its first publication in 1997, the Racah package fB] has grown considerably. In its early 
days, emphasis was put on simplifying those Racah expressions which included summations 
only over the Wigner n — j symbols [cf. Ref. 0, Figure 1]. Apart from such algebraic manip- 
ulations of expressions, however, recent developments also concerned an improved numerical 
support of a few standard tasks from the theory of angular momentum. Following this line, 
a number of basic quantities for evaluating (many-electron) matrix elements have been pro- 
vided recently 7 and will be further enhanced in the future. With the present extension to 
the Racah program, we now facilitate also the access to LS—jj transformation matrices and 
to the transformation of general atomic and configuration symmetry functions at different 
level of complexity. 

Obviously, any automatic transformation of atomic states must enable the user with a quick 
and simple access to the underlying symmetry functions. To support the definition and 
communication of such functions, two basic terms from the atomic shell model play a central 
role: atomic shell states and their successive coupling which, finally, leads to the definition of 
configuration state functions ( CSF). In order to simplify the handling of these functions, four 
auxiliarity procedures have been designed and introduced into the Racah program with the 
intention to keep all necessary information about a subhsell state or a CSF close together. 
These procedures such as csf_jj(), shell_jj(), . . . are defined separately for each coupling 
scheme, i.e. for LS— or jj— coupling, and, hence, could be easily extended to include other 
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coupling schemes in the future. The procedure Racah_set_coupling() also differs from most 
other commands in that it just 'assigns' the coupling order LS or SL of the orbital angular 
momenta and the spins to a global variable; this particular procedure must therefore be 
invoked prior to any other command which deals with LS — jj transformations. 

Although a simple notation has been worked out in order to construct CSF of any complexity, 
an automatic transformation of atomic states is currently supported only for configuration 
states with up to two open shells in LS— and up to four open subshells in jj— coupling, 
respectively. This limitation has arose from the number of recoupling coefficients which 
occur in any transformation and which grows rapidly if more open shells get involved. In 
the future, the program might be extended to more complex shell structures if demands arise 
from our work or from the side of users. 

At user's level, only a very few procedures have to be known in order to obtain the LS — jj 
transformation matrices or to transform some atomic or configuration states. Table 1 displays 
a short list of the main commands. Owing to the rapid increase of the complexity of many 
expressions, however, a (much) larger number of procedures had to be implemented at a 
lower level of the program. Many of them make use of previous development such as the 
routines for obtaining reduced coefficients of fractional parentage or for evaluating recoupling 
matrices. Since all transformation coefficients are evaluated directly to their numerical (i.e. 
either algebraic or floating-point) values, no additional data structures needed to be defined 
for the present work. 

The explicit expression (|14|) of a transformation matrix for two open shells showed that these 
coefficients are always reduced to the transformation of individual subshell states 



These matrix elements are stored internally in the program using the format 

[Q, L, S, J_big, N_l, Q_l, J_big_l, Q_2, J_big_2, factor, nom, denom] 
for 1 = 1, 2 and 

[w, Q, L, S, J_big, N_l , Q_l, J_big_l, Q_2, J_big_2, factor, nom, denom] 

for I = 3, as their recursive computation (from the matrix elements of the corresponding 
parents states) was found too slow for practical purposes. In this representation, the value 
of each LS — jj transformation matrix elements is given by factor x J™^ m ; they are kept 
for all occupation numbers N = 1, 2/ + 1. 

With the increased size of the Racah program, a new program structure became necessary. 
Following recent suggestions by the Maple standard, therefore, the program is now divided 
into the two modules Racah and Jucys which must be loaded separately by the withQ 
feature of Maple. While the Racah module now contains all procedures for the set-up 
and manipulation of Racah expression IHJ , the standard quantities from Ref. and the 
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present implementation of the LS — jj transformation matrices are incorporated into the 
module Jucys. Of course, the use of modules also helps to keep all low-level procedures 
invisible to the user. The Racah package is distributed as a tar (-xvf) file Racah2002 . tar 
which, apart of the source code and module libraries in different Maple versions, includes 
a Read. me for the installation of the program as well as the document Racah-commands.ps. 
This document provides the definition of all data structures of the Racah program as well 
as an alphabetic list of all user relevant commands. The code can be down-loaded also 
from our home page via the world-wide-web (http://www.physik.uni-kassel/fritzsche). For 
most commands, moreover, there are on-line help pages available which are distributed and 
maintained together with the code. 



4 Examples 

A few examples from atomic shell theory are displayed below to illustrate the application of 
the LS—jj transformation matrices from section 2. Beside of the computation of a particular 
matrix elements, we briefly explain how, for instance, the user can generate a (full) LS — jj 
transformation matrix for a half-filled /—shell, i.e. the subspace which is spanned by the 
|/ 7 wuLS) subshell states. We also show explicitly how, for C 2+ ions, the two low-lying 
ls 2 2s2p J = 1 levels from a (jj— coupled) multiconfiguration Dirac-Fock calculation can be 
easily transformed into a more appropriate LS— coupled basis. 

Let us start with the (numerical) evaluation of the LS — jj transformation matrix element 2 



(/ 3 (w = l,u = 3) 2 J 15/2 | f 7 % (t= 3) 15/2) 



for a partially filled /—shell, i.e. for the two (sub-) shell states with quantum numbers I = 3, 
N = 3, w = 1, v = 3, L = 7, S = 1/2, and + j= 7/2, JV= 3, v= 3, J= 15/2, respectively. 
Since N= in this case, no / 5 ^ 2 vJ^j subshell state occurs in the notation. For these 
quantum numbers, we obtain the transformation coefficient simply by typing 

> Racah_set_coupling(LS) ; 

> T := Racah_matrix_LS_j j (shell_LS(3 , 3 , 1 , 3 ,7 , 1/2) , 

shell_jj (-4,3,3, 15/2) ) ; 

T := .586845597 . 



As previously, the same result can be obtained also in algebraic or prime-number repre- 
sentations if the proper keywords algebraic or prime are added to the parameter list. In 
prime-number notation, for instance, the transformation coefficient 

> T := Racah_matrix_LS_j j (shell_LS(3 , 3 , 1 , 3 ,7 , 1/2) , 

shell_j j (-4,3,3, 15/2) ,prime) ; 

2 The capital J on the left-hand side denotes the spectroscopic notation for the orbital angular momentum 
L = 7 and should not be confused with the total angular momentum of this shell state. 
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T := [1, -3, 3, 1, -2] 

is returned as a list of (the first few non-zero) integer powers 

[do, d\, 0,2, as, «4, «5> «6> a 7, as, ag, aio, an] • 

of the prime numbers p\ = 2, p2 = 3, P3 = 5, p^ = 7, p$ = 11, ^6 = 13, p-j = 17, = 19, 
pg = 23, pio = 29, pn = 31 from which the actual values is obtained as 

«o ( Tip? ) 2 • (17) 



That is, the result [1, -3, 3, 1, -2] is just equivalent to the value ^ y ~ 0.586845597. 
Several tabulations have been published over the years |5l ll6llT7] to the LS—jj transformation 

/ M -N_- +N + + \ 

matrices ( I aLS J (k uJ, k vJ)J ) for different — partially filled — shells and using 

different phase conventions. With the present extension to the Racah package, we now 
provide a much simpler access to these transformation matrices which can be adopted to the 
actual requirements of the user. As a second example, therefore, we display how one can 
create an 'electronic table' (for any shell with I < 3) within only a few lines of Maple code. 
Figure 1 shows the necessary code for a half-filled / 7 shell which can easily be modified and 
extended for other shells as well. The printout from this example is shown below in the Test 
Run Output. 

We now extent our examples to the transformation of configuration state functions or even 
atomic states as they frequently appear in standard (relativistic) computations. To deal with 
a simple case, let us consider the two lowest ls 2 2s2p J = 1 levels of C 2+ ions which attracted 
a lot of recent interest in the diagnostics of stellar atmospheres [15] , In a single-configuration 
approximation, these two levels are written in terms of only two jj— coupled CSF 

|* a ) = ai(a) |7i J = 1) + a 2 (a) | 72 J = 1) (18) 

where a = 1,2, and the configuration states 

| 7 iJ = l) = \(2s\ /2 (y = l) 1/2; 2p\ /2 (i/ = 1) 1/2) J = l) 

172 J = 1) = | (2*1/2 iy = 1) 1/2; 2^ /2 {y = 1)3/2) J = l) ' 1<M 

are derived from the coupling of the two valence electrons in an open s— and p— shell, respec- 
tively. In this notation, we omit the Is 2 core since it does not take part in the coupling of the 
shells or in according transformations. The mixing coefficients {aj (a)} in Eq. ()18j) can be 
obtained from either a multiconfiguration Dirac-Fock (MCDF) or configuration interaction 
calculation. By using, for instance, the well-known Grasp92 code jlUj . we find for the two 
J = 1 levels the expansions 



l*i> 

1*2) 



= 0.8170 |7i J = 1) + 0.5767 \ j 2 J = 1) 

= -0.5767 I71 J = 1) + 0.8170 (72 J = 1) . (20) 
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Racah_set_coupling(LS) ; Racah_set_coupling_scheme (LS_quasispin) ; 

1 := 3; N := 7; jl := 5/2; j2 := 7/2; tl := Racah_subshell_term_LS(l , Q_int) ; 

kappa_l := 3; kappa_2 := -4; Racah_set_coupling_scheme ( j j_quasispin) ; 

for J from 1/2 to (N*l)+l/2 do lprint ( " J=" , J) ; 

for i from 1 to nops(tl) do 

if abs(tl[i] [4]-tl[i] [5])<=J and J<=tl [i] [4] -tl [i] [5] and type( J+tl [i] [4] +tl [i] [5] , integer) then 
si := shell_LS(l, N, tl[i][2], 2*l+l-2*tl [i] [3] , tl[i][4], tl[i][5], check); 
lprint (Racah_shell_print (si) ) ; 
for Nl from to N do 

Ml := (Nl-(2*jl+l)/2)/2; N2 := N-Nl; M2 := (N2-(2*j2+l)/2)/2; 

if type (Ml, integer) then 

t2 := Racah_subshell_term_j j (j 1 ,Q_int) ; t3 := Racah_subshell_term_jj (j2,Q_int) ; 

elif not type (Ml, integer) then 

t2 := Racah_subshell_term_j j (j 1 ,Q_half int) ; t3 := Racah_subshell_term_j j (j2,Q_half int) ; 
f i; 

for i_l from 1 to nops(t2) do 
if abs(Ml) <= t2[i_l][3] then 

s2 := shell_jj (kappa_l, Nl, (2*j 1+1) /2-2*t2 [i_l] [3] , t2 [i_l] [4] ) ; 
for i_2 from 1 to nops(t3) do 

if abs(M2) <= t3[i_2] [3] then 
if abs(t2[i_l] [4] - t3[i_2][4]) <= J and J <= t2[i_l][4] + t3[i_2][4] then 
s3 := shell_jj (kappa_2, N2, (2*j2+l) /2-2*t3 [i_2] [3] , t3[i_2][4], check); 
result := Racah_matrix_LS_j j (si, s2,s3,J, prime) ; 
lprint (Racah_shell_print (s2) ,Racah_shell_print(s3) .result) ; 
end if ; end if ; 
end do; end if; end do; end do; end if; end do; end do; 

Figure 1: Generation of the LS — jj transformation matrix for a half-filled / 7 shell. The 
beginning of this table is shown in the Test Run Output below. 

which clearly illustrate that a (pure) jj— coupling scheme is inappropriate for the present 
example. 

A (much) more appropriate representation might be obtained in LS— coupling. However, 
before we transform the two atomic states 1^1,2) into such a representation, we first show 
the transformation of a single configuration state, say, I71 J = 1) . The quantum numbers 
of this state can be read off directly from its definition in (|19|) . In the Racah program, we 
may enter this CSF as 

> CSF_1 := csf _j j (shell_j j (-1,1,1, 1/2), shell_jj (1,1,1, 1/2), 1, check); 

CSF_1 = csf_jj(shell_jj (-1,1,1, 1/2), shell.jj (1,1,1, 1/2), 1) 

where use is made of the two auxiliarity procedures csf_j j () and shell_j j (). These proce- 
dures return the input (basically) unevaluated but help facilitate the communication with and 
within the program (see Appendix B for further details about these commands). By defin- 
ing first (again) the coupling sequence for the LS— subshell states, we obtain the (complete) 
expansion of the CSF I71 J = 1) by 

> Racah_set_coupling(LS') ; 

> Racah_transform_csf("j j->LS" ,CSF_1, print) : 

".5773502693 * | (s~l nu=l, ~2S; p~l nu=l, ~2P) ~2P_1 >" 
".8164965809 * | (s~l nu=l, ~2S; p~l nu=l, ~2P) ~3P_1 >" 
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where, in spectroscopic notation, the first lines represents the l P\ (L = 1, 5 = 0, J = 1) and 
the second the 3 Pi component. In the last input line, the keyword print causes the procedure 
to 'print' the result (and to return a null expression) while, otherwise, the same result is 
returned in a list structure [csfiQ, ci, csf2(), C2, ...], as suitable for further manipulations. 

Having a simple access to the transformation of configuration states, we are now prepared 
to transform also the (full) atomic states. For the wave functions l^i^) of the two J = 1 
levels, this is simply achieved (as before) by assigning the CSF I72 J = 1) from ()19|) also to 
some variable CSF_2 and by carrying out the transformation explicitly: 

> CSF_2 := csf.jj (shelljj (-1,1,1,1/2) ,shell_jj (-2,1,1,3/2) ,1) ; 

> Racah.transf orm_asf ( " j j ->LS " , CSF_1 , . 8170 , CSF_2 , . 5767 , print ) : 

".9426 * |(s~l nu=l, ~2S; p~l nu=l, ~2P) ~1P_1 >" 
".3341 * I (s~l nu=l, ~2S; p~l nu=l, ~2P) ~3P_1 >" 

and 

> Racah.transf orm_asf ( " j j ->LS " , CSF_1 , -0 . 5767 , CSF_2 ,0.8170, print ) : 

" .3341 * I (s~l nu=l, ~2S; p~l nu=l, ~2P) ~1P_1 >" 
"-.9426 * |(s~l nu=l, ~2S; p"l nu=l, ~2P) ~3P_1 >" . 

Apparently, while the ground state level l^i) is a 1 P\ level, ^2) represents the 3 Pi level. 
Again, the print flag is used and ensures that the results are printed to screen and are not 
returned in terms of a list structure. 

The command Racah_transform_asf() provides a very flexible access to the transformation of 
atomic states; any number of CSF along with their corresponding mixing coefficients o« can 
appear in the parameter list. Moreover, the same syntax applies for this procedure if an 
atomic state is given in LS— coupling, 

|*«) = E < a ) \^ LSJ ) > ( 21 ) 

i 

and should be transformed into jj— coupling. With just two minor differences: The string 
"LS->jj" has to be used, instead, and the CSFj in the parameter list must represent proper 
LS*— coupled configuration state functions. Although the program is currently limited to two 
open shells, of course, the same syntax could be used for more complex shell structures or if 
other coupling schemes are to be incorporated into the program. Moreover, the construction 
of the CSF from the successive coupling of subshell states will help tackle more enhanced 
tasks in the future, such as the computation of angular coefficients for (non-scalar) tensorial 
operators of rank K. 
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5 Evaluation of transformation matrices 



We now return to the evaluation of the transformation coefficients IjllJI from section 2 which 
we could easily write down in this bra-ket notation. However, to carry out any transformation 
of configuration or atomic states explicitly, these matrix elements must be simplified to a 
computationally suitable form. This is achieved by the recoupling of the angular momenta 
which, in a number of steps, enables us to bring them into an equal sequence on the left- and 
right-hand side of the transformation matrix Therefore, any transformation matrix can 

always be expressed in terms of several recoupling coefficients and an appropriate number of 
transformation matrices (J7J), i.e. 

/ N , I-?*-- + + \ \ 

( I aLS J | \ k uJ, k vJ \ J ) , 

where a single matrix occurs for each open shell in the construction of the symmetry-adapted 
functions. 

For LS — jj transformations and the standard order (|TU)) of the jj— coupled subshell states, 
two steps are required for the recoupling of the angular momenta. The first step (i) arise 
from the recoupling of the total subshell orbital angular momenta Li and spins Si 

(((((ii, L 2 )L 12 , L 3 )L 123 , ...)L, (((Si, S 2 )Si 2 , S 3 )S 123 , ...)S) J \ 

((((Li, Si)Ji, (L 2 , S 2 ) J 2 ) J 12 , (L 3 , S 3 ) J 3 ) J123-) J) , (22) 

in order to obtain the total subshell angular momenta Jj. In the second step, then, (ii) these 
total angular momenta are brought into their standard order (jlOj) by 

((((Jl, Jl)Jl, (J 2 , J 2 )J 2 )Jl 2 , (J 3 , J 3 )J 3 )Jl 23 , ...) J I 

((((Ji, Ji) Ji, J 2 , )J{ 2 , J 2 ) J 12 , J 3 )J{ 23 , J 3 ) J 123 , ... ) j) . (23) 



Further steps in the recoupling of angular momenta may arise if the jj— coupled CSF are not 
defined in standard order (|lUj) or if more elaborate coupling schemes occur. 

Recoupling coefficients of type (|22j) and (|23|) can be evaluated by means of the Racah 
program [S]. In the following, we demonstrate this recent progress with the simplification of 
these coefficients for configuration states with two open shells in LS— coupling 
(l^aiLiSi, I2 ' 2 ' a.2L2S2)LS J^j. In this case, the first recoupling coefficient (f2*2l simplifies to 
{((Li, L2)L, (S\, S2)S)J\ ((Li,Si)Ji, (L2, S2) J2) J) and is evaluated interactively by 

> rcc_l := Racah_set(recoupling('<((Ll,L2)L, (S1,S2)S) J| 

((LI, SI) Jl, (L2,S2) J2) J>')) : 

> rcc_l := Racah_evaluate (rcc_l) : 

> Racah_print (rcc_l) : 
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— > 



(-2 LI + 2 S2 + 2 SI + 2 L2 + 2 J) 



(-1) 



V2J2 + 1V2J1 + 1V2S + 1V2L + 1 
w9j (L2.L1.L, J2, Jl, J,S2,S1,S) 



which can be re-written as 



<((Li, L 2 )L, (S x , S 2 )S)J | ((Li, Sj) J x , (L 2 , S 2 ) J 2 ) J) 

Li Si Ji 



[Ji,J 2 ,J,S] { 



J 2 S 2 J 2 
L S J 



(24) 



by using the symmetries of the Wigner 9 — j symbols and the fact that 2J + 2 mJ_ can be 
added without any change in the overall phase of the expression. Using similar lines, we can 
evaluate the second coefficient ()23[) for the given case of four open {jj— coupled) subshells 
{((Jmi,J P i) Ji, (Jm2,J P 2)J2) J\ ((( J ml , Jpl) Ji, Jm2)J'i2, Jpi) J ) in expression JT3J) 

> rcc_2 := Racah.set (recoupling( ' < ( ( Jml , Jpl) Jl, (Jm2, Jp2) J2) J| 

( ( ( Jml , Jpl) Jl , Jm2) J12_p , Jp2) J> ' ) ) : 

> rcc_2 := Racah_evaluate (rcc_2) : 

> Racah_print (rcc_2) : 



(-Jp2 + 2 mJ_ - Jl + 2 Jml + 2 Jpl - J - Jm2 + 2 J2) 



(-1) 



^2J12_p + 1V2J2+ 1 

triangle ( Jl , Jml , Jpl) 
w6j ( J , Jl , J2 , Jm2 , Jp2 , J12_p) 



which leads to the result 



{{{Jml, Jpl)Jl, {Jm2, Jp2)J2)J \ {{{Jml, JpljJl, Jm2)J\2, Jp2)J) 



(-i) Jm2+jp2+Ji+j v^y { / 

L Jp2 



J Jm2 Ji 



J 



12 
J 2 



(25) 



When we combine the two expression (1241) and (|25j). we arrive at the expression (|14[) for the 
transformation matrix of configuration states with two (open) shells in LS— coupling; their 
complete expansion in terms of a jj— coupled basis can be written 



(Z^aiLiSi, l^ 2 a 2 L 2 S 2 )LSJ 



Na, 



E 

,--++--++ 

Ni_ N 2 _ Ji JiJ xri v \J \viJ \viJ iviJ i 



y/[Ji,J2,L,S\ 



L\ Si Ji 
L 2 S 2 J 2 
L S J 
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l_ - + (Ni-N ± ) + + 



x ( (k 1 v\Ji, k x v\J\)J\ l^ 1 aiLiSi J\ 



x I (k 2 v 2 J 2 , k 2 v 2 J 2 )J\ l 2 2 a 2 L 2 S 2 J 2 



+(N 2 -N 2 ) ± + ^ ^ 
_N t _ - ,iVi, + _N 2 _ - ,Na 



+ - + 



(((«! V\J\, K x V\J\)J\, K 2 v 2 J 2 )J[ 2 , n 2 v 2 J 2 )Jj (26) 

or vice versa, if a summation is carried out over all intermediate angular momenta in 
LS— coupling. 

6 Summary and outlook 

A set of additional commands to the Racah program now facilitates the transformation 
of symmetry-adapted functions with quite different complexity from jj— to LS— coupling 
and vice versa. For such transformations, all partially filled (sub-) shells with / < 3 (i.e. 
up to /—electrons) are supported and, hence, the program extends the previously available 
tabulations and implementations considerably. In the study of atomic spectra, for example, 
the new version of the Racah program may help identify atomic and ionic levels as obtained 
from relativistic calculations in jj— coupling. Apart from the analysis of the valence-shell 
spectra, a reliable classification of the level structure is crucial, in particular, for the study of 
inner-shell processes, where the creation of additional vacancies often gives rise to a large (or 
even huge) number of possible states; in practise, however, only a very few levels are typically 
involved in some process but need first to be recognized, of course. The present extension to 
the Racah program can help to implement such transformations also directly into available 
atomic code, a project which is currently under work for the Ratip package jllj . 

With the implementation of rather abstract data structure such as shell and configuration 
states [cf. the auxiliarity procedures shell_LS() , csf_LS() , . . .], we also provide a pow- 
erful notation for more advanced tasks. For open-shell atoms and ions, for instance, a long- 
standing problem concerns the computation of the angular coefficients for effective n— particle 
operators as they occur in many-body perturbation theory. Here, the given notation for shell 
and configuration states can help decompose general matrix elements automatically. Another 
task concerns the optimal classification of atomic levels to assist the interpretation of atomic 
data and to improve the data base on energy levels and transition probabilities for the large 
user community of atomic data. For this aim, further coupling schemes need to be imple- 
mented in the future. 

Acknowledgement: We like to thank T. Inghoff for valuable suggestions and for help with 
the installation of the code. 
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Appendix A: Recurrence relations for the subshell transforma- 
tion coefficients 



The transformation coefficients ( I aLS J 



+ N + + \ 

(k vJ, k vJ)J ) between the (sub-)shell 



states from different coupling schemes are the building blocks for all transformations. For a 
given single-shell configuration of N equivalent electrons (I ), these matrices can be expressed 
recursively in terms of the transformation matrices for N — l equivalent electrons, i.e. in terms 
of the transformation of the corresponding parent states. The recurrence relations for the 
subshell transformation coefficients therefore include the coefficients of fractional parentage 



l N aLS 



x 1 (a'L'S') I ) and a proper recoupling of the angular momenta 



l N aLS J | (k vJ, k uJ)J 



_n_- _^N_ 



= V\L,S}/N ( lN aLSWl"- 1 (a'L'S') l) £ [J'j 



a'L'S' 



V I L 
N\j-,J]{S' s S 
J' 3- J 



V J 



x 3^ (v J ) , 



j N - vJ ) ( I"' 1 a'L'S' J'\(k vJ,kvJ )J' 



_{N-1)_I-' ,N,+ 



X [ 3+ (v J)j 



V I 






S' s 


:! 




J' 3, 




V J 


'+' 


+ 


+ \ / 


J)3+ 





+' 



J \j' J J' 

\j 3+ J 



J ) ( l( N -V a'L'S' J' | (ifuJ, t N J' 



I aLS J | (kpJ, kvJ) J ) = 



1 



(l + (-l) L+S ) y/\j-,j+,L,S\ 



They can be applied to the transformation of any subshell state by starting from 

I I L 

s s S 

. 3- 3+ J 
and 

2 ( I I L 

l 2 aLSJ\ttj\ = J(l + (-l) I+s )(l + H) J )[j±]/M S s S 

[j± 3± J 



(27) 



(28) 



(29) 



In the present work, the recurrence relations (|27|) have been utilized to generate the trans- 
formation matrices for all partially filled shells with I < 3, i.e. up to /—electrons, and for 
occupation numbers N = 1, 2, 21 + 1. In this process, we made use of the coefficients 
of fractional parentage (l N aLS l N ^ (a'L'S') /) in LSJ~ and (j*' 1 (y 1 ' J')j j N i/j) in 
j j— coupling, which were implemented earlier into the Racah package 1:,. To keep the (cur- 
rent) transformation of the symmetry-adapted functions feasible in time, these coefficients 
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are stored internally in the program [cf. section 4]. For all occupation numbers larger than 
21 + 1, i.e. N = 2Z + 2, 4/ + 1, we make use of the symmetry relation |18| 



l N avLSJ\ (k uJ, k vJ) J 



_N 



N 



A41+2-N a „LSJ\li 2i ~ +1 ~ m vj, ^ + >- m ti)j) . (30) 



which is easily derived from the following two properties of the coefficients of fractional 
parentage [cf. Eq. (15) in Ref. QU and Eq. (9) in Ref. [20]] 



(Vm-iv (a'u'L'S')l 



l 4l+2 ~ N avLS 



_l\S+S' '+L+V '-l-±+±(v+v' -1) 

(iV + 1) (2L' + 1) (25 , + l) 
(41 + 2 - N){2L + 1){2S + 1) 



l N (a'u'L'S')l U^'avLS) (31) 



JV+l 



and 



(f j ~ N (u ! J')3 \\j 2j+1 ~ N uJ) 



(jV + l)(2J' + l) 
(2j + l-N){2J + 



T) ) (j N {v'J')3 \\j N+1 uJ) . (32) 



For further details about the properties of the cfp coefficients and the subshell transformation 
matrices Q, see Refs. ^H^U] and Ref. (Sj, respectively. 



Appendix B: New commands for the Racah package 

The commands of the present extension to the Racah program can be described fairly 
independent from previous parts. Below, we briefly explain those procedures which have 
been added and which are of interest for an interactive use of the LS — jj transformation 
matrices. This provides a short description of the input and output of the procedures to 
facilitate also the understanding of our examples in sections 4 and 5; as previously, we follow 
the style of the former Maple Handbook ^1] . A more detailed description of all the presently 
available commands of the Racah package (at user's level) is distributed with the source 
code in the file Racah-commands .ps. 

As introduced earlier in the text, the terms of a subshell state and a configuration state 
function play a key role in the transformation of (coupled) states and the evaluation of 
matrix elements for open shells. They form the basic entities in dealing with such tasks 
and are often used to describe the input and output of (many) commands. To facilitate the 
handling of these 'atomic states' (i.e. the communication with and among the procedures of 
the Racah program), we first introduce a number of auxiliarity procedures for these coupled 
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states in LS— and jj— coupling. Although several tests are made on the particular input 
of these procedures, they basically return unevaluated and, thus, serve mainly for keeping 
necessary information together. As seen from their names and list of parameters, these 
auxiliarity procedures are designed so that further coupling schemes can be easily added 
later as according requirements arise. Since these procedures frequently occur during input 
and output and no actual manipulation is made, also no prefix Racah_ has been added. 

B.l Auxiliarity procedures 

Presently, procedures for the notation of shell-states and (atomic) CSF are provided in jj— 
and LS"— coupling. 

• csf_jj (shelLjj) 

Auxiliarity procedure to represent a configuration state function which is built from a 
single j j— coupled subshell state k n uJ^ or uk n v.J^ . 

Output: An unevaluated call to csf_jj (shelLjj) is returned. 

Argument options: (shelLjj, check) to check, in addition, that the given quantum num- 
bers in shelLjj () give rise to a valid jj— coupled subshell state with j = 1/2, . . . , 7/2. £ 
(shell_jji,shell_jj2,J) to represent a configuration state function (CSF) of two jj— coupled 



subshell states 



(k^ViJi, k 2 N2 v 2 J 2 )J^ or (niK^ViJi, n 2 K^ 2 u 2 J 2 )J^ 



& (shell_jji,shell_jj 2 ,Ji 2 ,shell_jj3,J) to represent a configuration state function (CSF) of 
three jj— coupled subshell states {{n^viJi, ^ 2 2 v 2 J 2 )J\ 2 ^ k^ 3 v 3 J 3 )J^ or 

((niK^ViJi, n 2 K 2 2 u 2 J 2 )Ji 2 , n 3 K^ 3 u 3 J 3 )J^ . 

Jit (shelLjji,shell_jj 2 ,Ji 2 ,shell_jj3,Ji 2 3,shell_jj4,J) to represent a configuration state func- 
tion (CSF) of four jj— coupled subshell states 
(((k^ViJi, k 2 N2 u 2 J 2 )J 12 , k^^J^Ji^, k^ 4 v 4 Ji)J^ or 

(((rai/^ViJi, n 2 K 2 2 u 2 J 2 )J 12 , n 3 K^ :i u 3 J 3 )J 123 , n±K^v± J 4 ) . 
Additional information: All (given) quantum numbers in the parameter list must 
evaluate to proper integers or half-integers. 

See also: csf_LS(), shell_jj(), Racah_csf_print(). 
• csf_LS(shell_LS,J) 

Auxiliarity procedure to represent a configuration state function which is built from a 



single LS*— coupled subshell state 



l N vLSJ) or nl N vLSJ) . 



Output: An unevaluated call to csf_LS(shell_LS,J) is returned. 

Argument options: (shelLLS, J, check) to check, in addition, that the given quan- 
tum numbers in shell_LS() and J give rise to a valid LSJ— coupled subshell state 
with I = 0, . . . , 3. Jit (shell_LSi,shell_LS 2 ,L,S,J) to represent a configuration state 
function (CSF) of two LS'-coupled subshell states (/{^ViLiSi, l2 2 u 2 L 2 S 2 )LSJ^ or 

(m/fV^iSi, n 2 l 2 N2 ^L 2 S 2 )LSj) . 
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Additional information: All (given) quantum numbers in the parameter list must 
evaluate to proper integers or half-integers. 

See also: csf_jj(), shell_LS(), Racah_csf_print(). 
• shell_jj(kappa,N,nu,J) 

K, N vJ) for 



Auxiliarity procedure to represent a jj— coupled subshell state 
j = 1/2, ...,7/2. 

Output: An unevaluated call to shell_jj(kappa,N,nu,J) is returned. 

Argument options: ([n,kappa],N,nu,J) to represent a jj— coupled subshell state 
nK N v.J^ . Jft (kappa,N,nu, J, check) to check, in addition, that the given quantum 
numbers give rise to a valid jj— coupled subshell state; the program terminates with an 
proper error message if this is not the case. 

Additional information: All quantum numbers (except of n) must evaluate to proper 
integers or half-integers. & The relativistic angular momentum quantum number k = 
± (j + 1/2) for / = j ±1/2. H> The principal quantum number n is often not required for 
the transformation of subshell states but enters the notation, if different subshell states 
are coupled to each other or, in particular, in the evaluation of (most physical) matrix 
elements. X All occupation numbers must be in the range N = 0, . . . , (2j + 1). £ For 
N = 0, an (unphysical) subshell angular momentum j = —1/2 is formally allowed in 
order to facilitate the input for several procedures from Appendix B.2. 

See also: csf_jj(), shelLLSQ, Racah_shell_print(), Racah_tabulate(). 



shell_LS(l,N,nu,L,S) 



l N vLS) for 



Auxiliarity procedure to represent a LS— coupled subshell state 
I = 0, ...,2. 

Output: An unevaluated call to shell_LS(l,N,nu,L,S) is returned. 

Argument options: (l,N,w,nu,L,S) to represent a LS*— coupled subshell state 
l N wvLS^j for I = 3 and the additional quantum number w = 0, . . . , 10. 

Jit ([n,l],N,nu,L,S) or ([n,l],N,w,nu,L,S) to represent the LS— coupled subshell states 
nl N vLS^j or nl N wvLS^j , respectively. & (l,N,nu,L,S, check) to check, in addition, 

that the given quantum numbers give rise to a valid LS— coupled subshell state with 

I = 0, . . . , 2; the program terminates with an proper error message if this is not the 

case. 

Additional information: All quantum numbers (except of n) must evaluate to proper 
integers or half-integers. Jit The principal quantum number n is often not required for 
the transformation of subshell states but enters the notation, if different subshell states 
are coupled to each other or, in particular, in the evaluation of (most physical) matrix 
elements. Jit All occupation numbers must be in the range N = 0, . . . , 2(21 + 1). 

See also: csf_LS(), shelLjjQ, Racah_shell_print(), Racah_tabulate(). 
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B.2 Commands for LS — jj transformations 



This Appendix lists the commands for the transformation of coupled (subshell and configu- 
ration) states where we utilize the (auxiliarity) notation from the previous part. This enables 
us with a very compact but still flexible notation for the input and output of the individual 
procedures; for example, a notation like . . . , shell_LS a , shell_LSft , . . . means that the user 
may type explicitly . . . ,shell_LS(l a ,N a ,nu a ,L a ,S a ) ,shell_LS(l6,N5,nu ;) ,L fe ,S fe ) , . . . in 
the parameter list or first assign these (unevaluated) calls to shell_LS() to any variables, say 
wa, wb, and later only use these variables at input time: . . . ,wa,wb, . . . . To 'extract' the 
quantum numbers from these unevaluated calls, the command Racah_tabulate() is used. 

• Racah_csf_print(csf_jj) 

Returns a string of type " |...('ki ~Ni, nui, Ji'; 're 2 ~N2, nu2, J2') J12; •■•>" to facilitate 
the printout of jj— coupled CSF. The value of m is printed in spectroscopic notation 
such as d_3/2, L7/2, ... ; if, moreover, the principal quantum number n is given, a 
string like 3d_5/2"2, ... is returned. 

Output: A string is returned. 

Argument options: (csf_LS) to return " |...('li ~Ni, nui, "2Si + l, Li'; 1 2 ~N 2 , nu2, 
"2S2+I, L2') L12, S12; ...>". The values of 1, and Lj are printed in spectroscopic 
notation such as s, p, d, ... and S, P, D, respectively. For /—electrons {U = 3), the 
additional quantum number Wj is printed in parenthesis such as "|...(li"Ni, (wi) nui, 
A 2Si+l, Li; 1 2 ~N 2 , (w 2 ) nu 2 , A 2S 2 +1, L 2 ) L 12 , Si 2 ; ...>". 
See also: csf_LS(), csf_jj() and Racah_shell_print(). 



Racah_matrix_LS_jj(shell_LS,shell_jj ,shell_jj + ,J) 

Returns the LS — jj transformation matrix / l N vLSJ (re uJ, re vJ)J\ for all 



.N_- ,N,+ \ 



subshells with I = 0, . . . , 3 and the according k = I and re = — (I + 1). 



Output: A (floating-point) number is returned. 

Argument options: (shell_LS,shell_jj_,shell_jj + ,J, algebraic) to return the same element 
of the LS — jj transformation matrix but in algebraic form. 

& (shelLLS,shell_jj_,shell_jj + ,J,pnme) to return the same element of the LS — jj trans- 
formation matrix but in prime-number representation. 

4k (shelLLSi,shelLLS2,L,S,shell_jji_,shelLjji + ,Ji,shelLjj2-,Ji2,shelLjj2+,J) to return 
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the LS — jj transformation matrix 



JV, 



+ N 1+ + 



,N 2 _ 



+ N 2+ + 



^(/j Viii^i, l 2 2 v 2 L 2 S 2 )LSJ (((«! viJi, n l v 1 Ji)J 1 , k 2 v 2 J 2 )J l2l k 2 v 2 J 2 )J j 

Additional information: The subshell angular momenta and occupation numbers are 

not independent of each other; they must fulfill the relation / = j_ + 1/2 = j + — 1/2 

+ 

and N = N + N; the program terminates with an proper error message if this is not 

the case. & If the principal quantum number n is given, it must be the same for all (sub- 

) shells. X For two and more coupled subshell states, these relations and condition must 

hold for each group of subshell states on the lhs and rhs of the transformation matrix. 

£ For I = follows N = and k = 0; such (unphysical) subshell states are formally 

allowed in the Racah program but can also be omitted from the list of parameters 

above. £ The subshell shell_jj_ with k or shelLjj+ with k can be omitted from the 

+ 

list of parameters above if N= or N= 0. X For details about the prime-number 
representation, see Racah_calculate_prime(). 

See also: Racah_set_coupling(). 

• Racah_set_coupling(i5) 

Defines the use of LS*— coupling for the l N uLS^ subshell states. 
Output: A NULL expression is returned. 

Argument options: (SL) to defines the use of SL— coupling for the l N vSL^ subshell 
states. 

Additional information: The information about the current coupling scheme for the 
LS subshell states is kept in the global variable Racah_save_coupling_LS; its default value 
is Racah_save_coupling_LS=LS. & If the coupling scheme of the subshell states is to be 
changes, this procedure must be called before any transformation is made. 

See also: Racah_set_coupling_scheme(). 

• Racah_shell_print(shell_jj) 

Returns a string " k, ~N, nu, J" to facilitate the printout of jj— coupled subshell states. 
The value of Ki is printed in spectroscopic notation such as d_3/2, L7/2, ... ; if, 
moreover, the principal quantum number n is given, a string like 3d_5/2"2, ... is 
returned. 

Output: A string is returned. 

Argument options: (shelLjj, state) to return "|/c, ~N, nu, J>". & (shelLLS) to return 
"FN, nu, "2S+1, L" or "FN, w, nu, "2S+1, L". The values of 1 and F are printed in 
spectroscopic notation such as s, p, d, ... and S, P, D, respectively. & (shelFFS, state) 
to return " |FN, nu, A 2S+1, F>" or "|FN, w, nu, '2S+1, L>". 

Additional information: These strings facilitate the line-mode printout of (coupled) 
subshell states and CSF. 
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See also: shell_LS(), shell_jj() and Racah_csf_print(). 

• Racah_tabulate(shell_jj) 

Return a table with all denned quantum numbers of a jj— coupled subshell state. 

Output: A table T with entries T [n] , T [kappa] , T [N] , T [nu] , and T [J] is returned. 

Argument options: (shelLLS) to return a table T with all defined quantum numbers of 
a LS-coupled shell state; it has the entries T [n] , T [1] , T [N] , T [w] , T [nu] , T [L] , 
and T [S] . 

Additional information: If some quantum numbers such as the orbital quantum num- 
ber / is not defined, fail is returned for the corresponding entry. 

See also: shell_jj(), shell_LS(). 

• Racah_transform_asf("jj >LS",csf_jji,ai,csf_jj2,a2,...) 

Expands an atomic state function, which is represented in a jj— coupled CSF basis 



i*«> = E | CSF r j ) - 



k 

into a basis of LS'-coupled CSF, i.e. 

I*a> = E|cSFf 5) )c^. 

i 

Output: A list [ [csf_LSi,ci], [csf_LS2,C2], ...] is returned where csLLSj describes a 
LS'-coupled CSF and C; the corresponding mixing coefficient in the expansion. 

Argument options: ("jj->LS",csf_jji,ai,csf_jj2,a2,...,a/fi r e6raic) to return the mixing co- 
efficients in algebraic form. X ("jj->LS",csf_jji,ai,csf_jj2,a2,...,prmi) to print the expan- 
sion in line mode. One line is printed per term c, * CSF^ LS ^\, and a null expression is 



returned in this case. £ ("LS->jj",csf_LSi,ai,csf_LS2,a2,...) to expand an atomic state 
function, which is represented in a LS— coupled CSF basis 



l*«> = E |cSFf 5) ) 4 LS \a) , 



k 

into a basis of jj— coupled CSF, i.e. 



|$ Q ) = ]T Icsf 



An) 



Additional information: The subshell states of all jj— coupled CSF must be provided 
in standard order, i.e. if both jj— subshells with j = Z zb 1/2 occur in a CSF, they must 
always couple like 

_A^- + N + + 

(k vJ,k vJ)J 
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if they represent the first two subshells, and 



_N_- +N + + 
((..., K vJ)J, K uJ)J ...) 

otherwise. & For any expansion into a jj— coupled CSF basis, the subshell states also 
appear in standard order in the output. $ If the principal quantum number(s) n are 
given, they are transfered properly to the output but must be the same for each group of 
subshell states in the expansion of the CSF. £ For the use of a SL— coupled CSF basis, 
see Racah_set_coupling(). & The subshell k or k can be omitted from the arguments of 

csLjj if N= or N= 0. 

See also: csf_jj(), csf_LS(), shell_jj(), shell_LS(). 
• Racah_transform_csf(" jj >LS" ,csf_jj) 

Expands a jj— coupled CSF into a basis of LS— coupled CSF 

CSF = |CSF. (L5) ) a. 



Output: A list [ [csf_LSi,ci], [csf_LS2,C2], ...] is returned where csLLSj describes a 
LS— coupled CSF and c« the corresponding mixing coefficient in the expansion. 

Argument options: (" jj~>LS" ,csf_jj, algebraic) to return the mixing coefficients in al- 
gebraic form. & ("jj— >LS",csf_jj,prmf) to print the expansion in line mode. One line 
is printed per term Cj * CSF^ LS ^, and a null expression is returned in this case. 
& ("LS->jj",csf_LS) to expand a LS— coupled CSF into a basis of jj— coupled CSF 



CSF (L5)\ = CSF wm c 



/ 



-Cm)\ 



/ 



Additional information: The subshell states of all jj— coupled CSF must be provided 
in standard order, i.e. if both jj— subshells with j = I ±1/2 occur in a CSF, they must 
always couple like 

^N_- + N + + 

(k vJ,k vJ)J 
if they represent the first two subshells, and 

_iV_- +N + + f 
((..., K uJ)J, K vJ)J ...) 

otherwise where k= I and k= — (I + 1). & For any expansion into a jj— coupled CSF 
basis, the subshell states also appear in standard order in the output. 
& If the principal quantum number(s) n are given, they are transfered properly to the 
output but must be the same for each group of subshell states in the expansion of 

the CSF. & For the use of a SL— coupled CSF basis, see Racah_set_coupling(). & The 

+ + 
subshell k or k can be omitted from the argument of the csLjj if N= or N= 0. 

See also: csf_jj(), csLLSQ, shelLjjQ, shelLLSQ. 
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• Racah_transform_csf_jj_LS(csf_jj) 



Expands a single jj— coupled CSF into a basis of LS— coupled CSF 



CSfC">) = |CSF. (L5) ) Ci . 



Output: A list [ [csf_LSi,Ci], [csf_LS2,C2], ...] is returned where csLLSj describes a 
LS— coupled CSF and C; the corresponding mixing coefficient in the expansion. 

Argument options: (csf_jj, algebraic) to return the mixing coefficients in algebraic form. 

See also: csf_jj(), csf_LS(), shell_jj(), shell_LS(), Racah_set_coupling(), 
Racah_transform_csf(). 

• Racah_transform_csf_LS_jj (csf_LS) 

Expands a single LS"— coupled CSF into a basis of jj— coupled CSF 



Output: A list [ [csf_jji,Ci], [csf_jj2,C2], ...] is returned where csLjjj describes a 
jj— coupled CSF and Cj the corresponding mixing coefficient in the expansion. 

Argument options: (csLLS, algebraic) to return the mixing coefficients in algebraic 
form. 

See also: csf_jj(), csf_LS(), shell_jj(), shell_LS(), Racah_set_coupling(), 
Racah_transform_csf(). 
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TEST RUN OUTPUT 

> Racah_LS_j j_calculate_table() ; 
"LS-jj transformation matrices for f subshells with occupation N 
"J=", 1/2 

"f-7, (w=0) nu=5, ~6F" 



"f 


_5/2~l 


nu=l , 


5/2", 


"f_7/2 


~6 nu=2, 2", 


[1, 2, 


1, 0, - 


3] 


"f 


_5/2~2 


nu=2, 


2", " 


f_7/2~5 


nu=3, 3/2", 


[-1, 4, 


1, -1, 


-4] 


"f 


_5/2~2 


nu=2, 


2", "f_7/2~5 


nu=3, 5/2", 


[-i, o, 


o, -l, 


-3, 1] 


"f 


_5/2~2 


nu=2, 


4", "f_7/2~5 


nu=l, 7/2", 


[1, 3, 


0, 1, - 


3] 


"f 


_5/2~2 


nu=2, 


4", " 


f_7/2~5 


nu=3, 9/2", 


[-i, o, 


0, o, 


-4, 1, 1] 


"f 


_5/2~3 


nu=l , 


5/2", 


"f_7/2 


"4 nu=2, 2", 


[0] 






"f 


_5/2'3 


nu=l , 


5/2", 


"f_7/2 


~4 nu=4, 2", 


[1, o, 


1, 0, - 


3, 1] 


"f 


_5/2'3 


nu=3, 


3/2", 


"f_7/2 


~4 nu=2, 2", 


[1, 1, 


1, 0, - 


4] 


"f 


_5/2~3 


nu=3, 


3/2", 


"f_7/2 


"4 nu=4, 2", 


[0] 






"f 


_5/2~3 


nu=3 , 


9/2", 


"f_7/2 


"4 nu=2, 4", 


[1, 2, 


0, 2, - 


4, 1] 


"f 


_5/2~3 


nu=3 , 


9/2", 


"f_7/2 


~4 nu=4, 4", 


[0] 






"f 


_5/2~3 


nu=3, 


9/2", 


"f_7/2 


~4 nu=4, 5", 


[0] 






"f 


_5/2~4 


nu=2, 


2", "f_7/2~3 


nu=3, 3/2", 


[-1, 4, 


1, -1, 


-4] 


"f 


_5/2~4 


nu=2, 


2", " 


f_7/2~3 


nu=3, 5/2", 


[-i, o, 


o, -l, 


-3, 1] 


"f 


_5/2~4 


nu=2, 


4", " 


f_7/2~3 


nu=l, 7/2", 


[-1, 3, 


0, 1, 


-3] 


"f 


_5/2~4 


nu=2, 


4", " 


f_7/2~3 


nu=3, 9/2", 


[-i, o, 


0, o, 


-4, 1, 1] 


"f 


_5/2~5 


nu=l , 


5/2", 


"f_7/2 


~2 nu=2, 2", 


[-1, 2, 


1, 0, 


-3] 
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